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Abstract 

We present a convariant formulation for radiative transfer in curved 
space time and demonstrate some applications in the black-hole sys- 
tems. We calculate the emission from semi-transparent accretion tori 
around black holes, for opacity provided by the Fe Ka and K/3 lines 
and for opacity dominated by electron scattering. We also calculate 
the emission from radiative inefficient accretion flow in black holes 
with opacity provided by electron-positron annihilation lines. Finally 
we show shadows cast by accreting black holes with different spins and 
with different distribution of warm material around them. 



1 General Relativistic Radiative Transfer 

A black hole is characterised by a singularity and an event horizon. The 
presence of black holes in this universe is based on inferences, as no black 
hole has been observed directly. Now the question is: what does a black 
hole look like? Perhaps, it may be more meaningful to ask: what does an 
astrophysical black hole look like if we can image it, and can we see the 
shadow of its event horizon? 



Black holes are relativist ic objects. A convariant radiative transfer for- 
mulation is therefore needed to calculate the radiation from a black hole's 
vicinity. In the absence of scattering, the transfer equation reads 
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(Lindquist 1966; Baschek et al. 1997; Fuerst and Wu 2004, 2007; Wu and 
Fuerst 2008), where is the 4-momentum of the photon and is the 
4-velocity of the medium. The specific intensity X, absorption coefficient 
Xo and emission coefficient 770 are in Lorentz invariant forms. The left side 
of the equation describes the radiation propagation; the right side specifies 
how the radiation interacts with the medium. When the velocity fields of 
the emitters and absorbers are given, the transfer equation can be solved 
along the geodesies using a ray-tracing algorithm. 

Scattering deflects photons into the line-of-sight. In its presence the 
radiative transfer equation is an integro-differential equation. Fuerst (2005) 
used projection moment tensors to decompose the transfer equation into a 
set of moment equations. The first two moment equations are 
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respectively, where p is the density of the medium, is the absorption 
cross section, cTgc is the electron scattering cross section, <Sx is the source 
function, is the photon directional unit vector {n^ria = 1), £' = —k^u^^ 
— + — k^ /E^ and ^ = u^ria-fsn^ + n^Ua-fsu^- (Here, we assume 



that electron scattering is the prime scattering process.) 
specific intensity is 



The j-th order 



(4) 



(Fuerst 2005; Wu et al. 2006). The moment equations are not closed. How- 
ever, they would converge rapidly for certain choices of moment functions. 
In practical calculations, using a truncated set of low order moments is usu- 
ally sufficient. 






Figure 1: Surface brightness of a serai-transparent accretion torus around a Kerr 
black hole (a = 0.998), viewed at inclinations of 15°, 60° and 85° (left to right). 
The opacity is due to the Fe Ka and K/3 lines. The brightness is normalised, with 
all tori having equal maximum pixel brightness. 





Figure 2: Surface brightness of a semi-transparent electron scattering accretion 
torus around a Kerr black hole (a = 0.998), viewed at inclinations of 45° (left) and 
85° (right). Brightness normalisation is the same as that in Fig. 1. 



2 Imaging Black Holes 

Figure 1 and 2 show the images of accretion tori around Kerr black holes. 
The opacities are provided by the Fe Ka and Kf3 lines and by electron 
scattering respectively in the two cases. For the case in Figure 1, the non- 
scattering formulation (Eq. 1) is used in the calculations; for the case in 
Figure 2, the moment radiative transfer formulation (Eq. 2, 3 and 4) is 
used. The velocity field and density distribution in the tori are determined 
using a empirical velocity law gauged by magneto-rotational instability sim- 
ulations of accretion disks (see Fuerst and Wu 2004, 2007). The brightness 





Figure 3: Intensity images of emission lines from vertically Gaussian-distributed 
hot plasmas around a Schwarzschild black hole (top row) and a Kerr black hole 
with a = 0.998 (bottom row). The density distribution takes the form: n{r^ z) oc 



-3/2 



exp{—z /H ). The viewing inclinations are 75°, 45° and 15° (left to right 



columns). The opacity is provided by the electron-positron annihilation line. 



distributions of the tori in the images reflect the combined effects of gravi- 
tational and kinematic frequency shifts, Doppler boosting, and line-of-sight 
path length variations caused by gravitational lensing and frame dragging. 

Figure 3 shows the emission images of electron-positron annihilation lines 
from hot accretion flows. We adopt a vertically Gaussian density distribution 
n(r, z) (X z"^/^ exp(— z^/i/^), where the scale height H is determined by the 
thermal sound speed and the rotational Keplerian angular velocity (see Ball 
2008 for details). This density distribution is appropriate for hot radiative 
inefficient accretion flows. A similar density distribution was used in the 
black-hole shadowing calculations by Huang et al. (2007). The emission 
images of the Schwarzschild and Kerr black holes are distinguishable with 
the differences more obvious for larger viewing inclination angles. 

Figure 4 shows the shadows in a starry background sky cast by black 
holes with a Bondi-Hoyle accretion flows, i.e. the density distribution of 
material n(r) oc r~^/^. The emissions from the stars are gravitationally 



Figure 4: Shadows of a background sky cast by Schwarzschild black holes (left 
column), and Kerr black holes with a = 0.998 viewed at inclination of 90° (middle 
column) and 45° (right column). The optical depth is provided by a spherical inflow 
of warm accreting gas with a temperature T = 5 x 10^ K. The density distribution of 
the gas takes the form: n(r) ex r~^/^. The line-of-sight opacity r = for the black 
holes in the top row, and the mean line-of-sight opacity r = 5 for the black holes 
in the bottom row. The stars in the background sky are randomly generated and 
have the same brightness. The red line denotes the locations of geodesies crossing 
the event horizon. 

lensed, but the stellar emissions have no frequency shifts because the stars 
are at infinite distances. However, the stars appear to distribute unevenly 
in the images and the centroid of the shadow appears to shift when the 
black hole rotates. These phenomena happen for two reasons. The first 
is due to 'geometrical effect', which is caused by rotational frame dragging 
of the black hole. This can be seen easily when comparing the images in 
the top panels. The shadow is not circular and the centroid shifts if the 
black hole rotates. The distortion depends on both the black-hole spin 
and the viewing inclinations. The second is due to 'optical depth effect'. 
Gravitational lensing and frame dragging have different effects on the lines- 
of-sight for different impact parameters and for different photon incidence 
with respect to the black hole's spin axis (in the case of Kerr black holes). 



The integrated optical depths are not the same along different lines-of-sight. 
In the bottom panels of Figure 4, the emissions of the stars at the left side 
of the event horizon are more strongly absorbed than the emissions of the 
stars at the right side. 
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DISCUSSION 

THOMAS SCHWEIZER: Is it possible to measure the spin of a black 
hole and the last stable orbit by looking at the emission lines? 

KINWAH WU: In principle, yes. In practice, it is more complicated to 
deduce the black-hole parameters from the emission lines because of degen- 
eracy. The emission line profiles are not uniquely determined by the location 
of the last stable orbit of particles, and hence the black hole's spin. 



